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We prove that the necessary condition for a solid to be also a superfiuid is to have zero-point 
vacancies, or interstitial atoms, or both, as an integral part of the ground state. As a consequence, 
superfluidity is not possible in commensurate solids which break continuous translation symmetry. 
We discuss recent experiment by Kim and Chan [Nature, 427, 225 (2004)] in the context of this 
theorem, question its bulk supersolid interpretation, and offer an alternative explanation in terms 
of superfiuid helium interfaces. 



PACS numbers: 67.40.- 



67.80.-S, 05.30.-d 



Recent discovery by Kim and Chan 0,01 that solid ^He 
samples have a non-classical moment of inertia (NCMI) 
is a breakthrough result which has prompted renewed 
interest in the supersohd (SFS) state of matter. Early 
theoretical work by Andreev and Lifshitz 3] and Chester 
showed that solids may feature a Bose-Einstein con- 
densate of vacancies (or interstitial atoms) and thus pos- 
sess superfiuid (SF) properties. One may consider their 
work as establishing sufficient conditions for SFS. It was 
natural then to interpret mass decoupling in the torsion 
oscillator experiments as originating from small (about 
~ 1%) concentration of zero-point vacancies 0. 

However, the overwhelming bulk of previous experi- 
mental work (for review, see, e.g., 5]) indicates that va- 
cancies and interstitials in ^He crystals are activated and 
their concentration is negligible below 0.2 K . The most 
recent study @ looked at the density variations of solid 
^He between two capacitor plates and did not reveal any 
presence of vacancies. To deal with these facts an idea 
was put forward that strong exchange processes in quan- 
tum crystals may lead to superfluidity even in the ab- 
sence of zero-point defects 011L Mistakenly, this idea is 
attributed to Leggett's paper Leggett estabhshed a 
link between the SF response and the connectivity of the 
groundstate wavefunction and derived a rigorous formula 
for the upper bound on the superfiuid density, ps. Crys- 
tal defects and their relation to the connectivity was not 
discussed in Ref. 0. A separate issue is the "extremely 
tentative" order-of-magnitude estimate < 3 x 10~^ 
based on the exchange integral between '^He atoms Q, 
which, apparently, caused the misleading interpretation 
that interatomic exchange, on its own, may result in SF 
in the absence of vacancies ^3 • 

The central point of the discussion to follow is to an- 
swer the question whether superfluidity is possible in 
crystal structures with the number of atoms being com- 
mensurate with the number of lattice points and what 
are the necessary conditions for this to happen. 

Below we re-examine Leggett's work and show that 
it implies vacancies and/or interstitial atoms as a nec- 
essary condition for SFS in bosonic systems similar to 
'''He. Chester's "final speculation" that without them 
the solid state of matter is insulating P| proves to be a 



theorem. We present an alternative proof of the theorem 
using path-integral language in which the presence of va- 
cancies in the SF state is seen explicitly. It also provides 
a simple picture showing why exchange processes on their 
own do not lead to SF. Thirdly, we put forward an ar- 
gument based on the phase-particle-number uncertainty 
relation which relates SF, compressibility of pinned 
solids and vacancies (pinning is crucial to separate and 
suppress the contribution to the compressibility coming 
from the change of the lattice constant from that due to 
adding/removing particles to/from the bulk ^2)- The 
answer to the central question is then that SFS ground- 
states in commensurate solids have "zero measure" to be 
found in Nature, because they require an exact symme- 
try between zero-point vacancies and interstitials which 
is immediately broken by changing system parameters, 
e.g. pressure 13] . By excluding a bulk supersolid in- 
terpretation of the Kim & Chan results we are forced to 
look for an alternative explanation of their data based on 
the physics of disordered and frustrated ^He interfaces. 

As shown by Feynman the groundstate of the 

bosonic system has no zeros, g(xi,X2, ■ ■ ■ ,xn) 0. 
Moreover, in superfluids ^'g does not become macro- 
scopically small when one or several coordinates, say, 
xi,X2, ■ ■ - Xm, are taken around the system while other 
coordinates are kept fixed. This property (called connec- 
tivity by Leggett Q) is key for superfluidity, and is just 
another way of saying that topological off-diagonal long- 
range order (TODLRO) is required for SF 15, 16, IJ]- 
The requirement that connected be single-valued 
leads to the quantization of circulation and thus stabil- 
ity of persistent currents in samples with the cylindrical 
annulus geometry [^. 

To illustrate the point, consider a one-dimensional sys- 
tem of two identical bosons forming a bound (molecu- 
lar) groundstate, </?o(ki — ^2!), with localization length 
I. Naively, the first rotating state of the molecule on 
a ring of large circumference, L ^ ^, is written as a 
product of the the plane wave for the center of mass 
coordinate, R — {ri + r2)/2, times the bound state: 
ipi — e*^'^^/^ — f2\)- This expression, however, 

is not single-valued, because if ri or is taken around 
the ring we get ipi — > —(pi- The correct solution is to 
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replace ipo with ipo, which has a zero at |ri — « i/2, 
i.e. in the region where ipQ is exponentially small; at dis- 
tances |ri — -C i the two functions are almost identi- 
cal, \(po\ « (fQ. Note, that the energetic cost of creating 
a zero in this case is exponentially small and vanishes in 
the limit of infinite system size. 

The same considerations apply to the disconnected 
crystal state consisting of N bosons when decays to 
the macroscopically small value when any finite number 
of coordinates are taken away from their typical posi- 
tions in the crystal (other coordinates are kept fixed) and 
moved around the system. The first rotating state in the 
system with periodic boundary conditions can be written 
as 01 = e"^"^^^^^ ifo with R = X^i^i '''i/^ and (fo having 
zeros in regions where the modulus of ipo is exponentially 
suppressed and thus extra zeros do not cost finite (system 
size independent) energy. Clearly, the phase gradient cir- 
culation of (pi is ^ 1/-/V, and this system will not show 
the NCMI which is based on the impossibility of setting 
system in rotation with a macroscopically small veloc- 
ity or circulation. This should be compared with the 
first rotating state of the single-atom superfluid system, 
^(SF) _ gi27rX]fLi i-i/L with ipo > and phase gradi- 
ent circulation 27r. Now, making zeros in connected ipo is 
so costly energetically that the lowest energy state corre- 
sponds to the relatively high kinetic energy of rotation. 
For definiteness, we consider below only single-atomic su- 
perfiuids, but all considerations are readily generalized 
to the m-atomic case. We have essentially reproduced 
above the Leggett's argument that SF wavefunctions are 
necessarily connected. 

By definition, \'^Gi^iT^2, ■ ■ ■ ,xn)\'^ is the probability 
density to observe particles at the specified positions. We 
fix all coordinates except one, xi, and observe that in 
connected wavefunctions |^'g(xi)P remains finite when 
xi is taken arbitrary far from the initial position. For- 
mally, this property is identical to statistical properties 
of atomic configurations in classical crystals at finite tem- 
perature and was used by Chester to introduce vacancies 
in the ground state. This final correspondence was not 
elaborated in Ref. 0. 

How do we "visualize" vacancies/interstitials in the 
state of identical particles with strong exchange, espe- 
cially when the number of atoms coincides with the num- 
ber of lattice points? It appears that the common per- 
ception is that such solids do not have vacancies and 
interstitials, by definition, or else they are indistinguish- 
able from the standard zero-point fluctuations. Imagine 
a solid sample pinned by an (arbitrarily weak) external 
potential preventing it from moving as a whole. There 
is no problem in identifying lattice points using the av- 
erage particle density profile p{r) which is periodic in 
space. Consider now some typical spatial configuration 
of particle positions along with the lattice points and 
start the coarse-graining procedure of "erasing" the clos- 
est particle-lattice point pairs in the spirit of the spa- 
tial renormalization group treatment. As we progress 
towards mesoscopic length-scales all short-ranged zero- 



point fluctuations of atoms away from lattice points will 
be erased from the picture. The procedure continues un- 
til we have erased all pairs with sizes much smaller than L 
but much larger than all microscopic scales. We say that 
the crystal state has no vacancies and interstitials if the 
final coarse-grained configuration is empty. On the other 
hand, if the coarse-grained configuration still contains 
lattice points, or particles, or both, at arbitrary large 
distances, we say that it has crystal defects in it. The 
decimation procedure explains how vacancies are possi- 
ble in commensurate solids, and perfectly agrees with the 
conventional view of classical crystals at finite tempera- 
ture. For the commensurate solid with connected ^'g we 
may start with the perfect-lattice configuration of parti- 
cle coordinates and its empty coarse-grained picture, and 
then move xi arbitrary distance away to produce an im- 
age of the vacancy and interstitial. This will not result 
in the exponential suppression of the configuration prob- 
ability (in fact, such configurations will dominate in the 
normalization integral) . 

In the absence of exact interstitial/ vacancy symme- 
try the concentration of vacancies, Uy, and interstitials, 
TT-int, in the supersolid will be different, e.g. riy > n-mt, 
since broken continuous symmetry allows production of 
excess vacancies by making small changes in the lattice 
constant. [This mechanism is is not available in discrete 
models, and then Uy — rtjnt is possible.] In the transla- 
tion invariant system at T=0 one expects then ps = Ariy, 
just like in any other interacting bosonic system. 

Our second consideration is based on the path-integral 
formulation of quantum statistics 0| in terms of many- 
body trajectories, {xi(r)}, in imaginary time r e [0,/3] 
with periodic boundary conditions {xi{P)} = {xi{0)}. 
The most important superfluid characteristic of parti- 
cle trajectories, or world lines, is their winding num- 
bers, AI", a = 1,2, ... ,d. We assume periodic boundary 
conditions in space in all d-dimensions with linear sizes 
L" = L. To determine M" imagine a cross-section go- 
ing through point R perpendicular to the direction a 
and count how many times particles cross it from left to 
right, /c", and from right to left, fc". By definition, wind- 
ing numbers are M" = A:" — fc" . They are independent of 
the cross-section location R because trajectories are con- 
tinuous and periodic in imaginary time. The superfluid 
density is then given by 

p'^'^ = 2mTL^-'^ { M°'M'') , (1) 

where m is the particle mass. In d = 3, the superfluid 
density is finite in the thermodynamic limit, L — > oo, 
T — > 0, and T/L ^ 0, if the probability of having world 
lines with non-zero winding numbers in the ground state 
is close to unity. 

We now demonstrate that crystal states without zero- 
point vacancies are described by world-line configurations 
with M = 0, i.e. they are not superfluid. We start 
with the picture of a perfect crystal with particles tightly 
localized around equilibrium lattice points. In this state 
the trajectories are nearly straight lines with M = as 
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FIG. 1: Particle worldlines in different crystals at low temperature. The time axis is vertical. Dashed lines show the equilibrium 
lattice points, (a) Nearly classical crystal at low temperature; M = 0. (b) Insulating quantum crystal with large zero-point 
fluctuations and frequent particle exchange processes; M = 0. (c) Particle worldlines with non-zero winding number. 
I 



shown in Fig. 1(a). When tunnelling exchange processes 
are added into the picture the world lines are no longer in 
one-to-one correspondence with the lattice points. The 
nature of the exchange process, however, is such that 
when one particle leaves its equilibrium crystal point Ri 
and goes to point i?2, the other particle goes from R2 
to i?i (for pairwise exchange). Thus, the net current 
of world lines through any cross-section remains strictly 
zero. The same conclusion follows from the consideration 
of multiparticle exchange events 0, see Fig 1(b). 

Consider now a world-line configuration with non-zero 
winding number. Fig. 1(c). At any moment of imagi- 
nary time we consider the spatial configuration of par- 
ticle positions and lattice points and apply the coarse- 
graining procedure discussed above. Again, all short- 
lived and short-ranged exchange process and zero-point 
fluctuations will be erased once we pass several atomic 
distances, and in the insulating state the "movie" of the 
coarse-grained configuration evolution in time will show 
an empty "screen" for Figs. l(a,b). If M ^ 0, as in 
Fig. 1(c), there is no way for the renormalization proce- 
dure to erase all particles at all moments in time since 
topologically winding numbers correspond to particle tra- 
jectories moving continuously in the same spatial direc- 
tion and thus creating imbalance between particles and 
lattice points at arbitrary large distances. After all other 
particles are associated with lattice points and erased, 
the winding trajectory describes an interstitial-vacancy 
pair which separates over the distance of order L and 
eventually makes a closed loop around the system, see 
Fig. 13 For the statistics of such loops to give non- 
vanishing ( M"M'^) in the thermodynamic limit, they 
have to be typical in a given lattice structure. In other 
words, zero-point vacancies and interstitials are an inte- 
gral part of the groundstate. 

Our last consideration is based on the relation between 
the pinned compressibility, k, and zero-point vacancies. 
Compressibility can be calculated through the energies 
of states with one extra particle, Ejsi+i, and one extra 
vacancy, Em-i, as k = l/I^Ai?, where V is the system 
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FIG. 2; Evolution of the coarse-grained picture for the tra- 
jectory with non-zero winding number similar to Fig. 1(c). 
Filled and open circles show particle and lattice site positions 
correspondingly. Arrows indicate the direction of the particle 
number current. 



volume, and Ai? = -Bat+i — -Btv-i- Incompressible states 
have finite /S.E. On another hand, in a macroscopic sys- 
tem Ai? can be obtained by considering the energy in- 
crease by creating an interstitial-vacancy pair with ar- 
bitrary large separation between them (in pinned solids 
the notion of vacancy or interstitial is rigorously defined 
within the coarse-graining procedure as the presence or 
absence of perfect registry between particles and unit 
cells). Crystals without zero-point defects are gapped 
with respect to the interstitial-vacancy production (oth- 
erwise these defects would be an essential part of the 
groundstate) and thus are incompressible (if pinned) and 
vice versa. One step further, this implies that superfluid- 
ity and pinned compressibility come together and either 
one (including long-wave acoustic properties with addi- 
tional sound mode) can be used for the detection of the 
SFS state experimentally. This conclusion is in line with 
the famous uncertainty relation [Tll | between the phase 
of the superfluid order parameter, 0, and particle num- 
ber, A0AiV > 1/2. Because of this relation, one may not 
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introduce a well defined phase field for the incompress- 
ible state of matter which tends to completely suppress 
particle number fluctuations. 

We have little doubt that large activation energies 
for vacancies and intcrstitials in "^He measured down to 

1 K temperatures 5] will not radically change to near 
zero at lower temperatures, and that helium is a commen- 
surate solid at T = 0. Since it has no symmetry between 
the vacancies and interstitials, we conclude that there 
are no zero-point vacancies in bulk solid ^He. By ex- 
cluding superflow through the crystal bulk we are forced 
to look more closely at the superfluid properties of dis- 
ordered helium-substrate layers and frustrated interfaces 
between helium micro-crystallites. 

There are indications of strong disorder in experimen- 
tal system of Ref. The dependence of the super- 
fluid density on reduced temperature parameter t = 
{Tc — T) /Tc has little to do with the expected bulk super- 
fluidity t'^-^'^^ dependence. Instead, ps appears to vanish 
at Tc with zero derivative. Such a behavior can be mod- 
elled by a broad distribution of transition temperatures 
in the heterogeneous sample. This observation correlates 
with the gradual decrease of the decoupled mass with 
the increase of the torsion oscillator amplitude by orders 
of magnitude. Let us assume that a sample consists of 
micro-crystallites of linear size D with superfluid inter- 
faces of typical thickness d between them. The superlfuid 
fraction may be estimated then as Ps/ p ^ d/D. To have 
~ 1% of the superfluid mass coming from interfaces with 
d 10 A one will need crystallite sizes about a frac- 
tion of a micron. The variety of interfaces with different 



crystallographic indices provides a broad distribution of 
transition temperatures. 

One of the experimental mysteries is extreme sensi- 
tivity to the addition of ^He impurities at the level of 
JT-im ~ lOOppm. To minimize kinetic energy, light "^He 
atoms are likely to end up at frustrated interfaces, and 
then at the edges where different interfaces meet. This 
may increase ^He edge vs bulk concentration by a fac- 
tor as large as {D/dY and produce n^^^"'^ ~ 1. This 
will have a profound effect on the edge-connected inter- 
face superfluidity (edges then act as a disordered two- 
dimensional network of Josephson junctions). 

We are not aware of any systematic study what are the 
properties of interfaces between the '^He micro- crystals 
and helium ''He crystals on disordered substrates at el- 
evated pressures. Model simulations of domain wall 
boundaries in the checkerboard solid (obtained for in- 
teracting lattice bosons at half-integer filling factor) show 
that they remain superfluid deep into the bulk solid phase 
[20] . In the outlined picture three predictions are certain: 
(i) the superfluid fraction must strongly depend on the 
crystal growth process, (ii) the amount of ^He sufficient 
to suppress superfluidity scales Pj, and (iii) 

transition temperatures on interfaces do not depend on 
D and ps{T = 0). 
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